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Abstract 


We present two numerical approximations with non-uniform meshes to the 
Caputo—Fabrizio derivative of order a (0 < a < 1). First, the L1 formula is 
obtained by using the linear interpolation approximation for constructing 
the second-order approximation. Next, the quadratic interpolation ap- 
proximation is used for improving the accuracy in the temporal direction. 
Besides, we discretize the spatial derivative using the compact finite differ- 
ence scheme. The accuracy of the suggested schemes is not dependent on 
the fractional a. The coefficients and the truncation errors are carefully 
investigated for two schemes, separately. Three examples are carried out 
to support the convergence orders and show the efficiency of the suggested 
scheme. 
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1 Introduction 


In recent years, some approximations have been proposed for the fractional 
derivative of order a, such as Grtinwald—Letnikov, Lubich, and Caputo ap- 
proximations [15, 16, 24, 28, 30, 27]. The schemes that are proposed until 
now to discretize the Caputo fractional derivative have been limited to the 
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accuracy of order 2—a (0 < a < 1) on uniform meshes. One disadvantage of 
previously proposed methods is that when a = 1, its accuracy may lead to 
poor accuracy. Thus, in this point of view, the numerical solutions of high 
dimensional partial fractional differential equations require a large number of 
computations. Besides, from the truncation error estimate of the methods on 
uniform meshes verify that the accuracy is dependent on the fractional order 
a. Weare able to overcome these difficulties using Caputo—Fabrizio fractional 
derivative with a non-singular kernel. Afterward, we will obtain the second 
and third-orders accuracy in time that is independent of the fractional order 
Q. 

Let us consider the following time fractional diffusion and advection equa- 
tions, respectively: 


o? t 
eppu(x,t) = 982) + F2,0), reo, 0<t<T, 
0 t 
TDPu(X,t) = MEO 5 Hey, we, O<E<T, (1) 


in which (FD¢ is the ath Caputo—Fabrizio fractional derivative defined by 


t 
oD usb) = Ma) | u'(X, 8) exp (— a —)as, O<a<1, (2) 
where M(q) is a normalization function such that M(0) = M(1) = 1. 

In 2015, Caputo and Fabrizio [7] suggested a new definition of fractional 
derivative based on the exponential kernel. They considered two different 
representations for the temporal and the spatial variables. It is important 
and interesting that this approach describes the behaviour of classical vis- 
coelastic materials, electromagnetic systems, thermal media, and so on. An- 
other interesting property of this definition is that it opens up new avenues 
in the mechanical phenomena, related to plasticity, fatigue, damage, and 
electromagnetic hysteresis [7]. 

Recently, studies of Caputo—Fabrizio fractional derivatives have been car- 
ried out by some authors. Authors of [6] investigated the existence of a 
solution for two high-order fractional integro-differential equations including 
the Caputo—Fabrizio derivative. Atangana and Alqahtani [4] considered a 
numerical approximation of the space and time Caputo—Fabrizio fractional 
derivative in connection with ground water pollution equation. In [18], the 
authors presented a Crank—Nicolson finite difference scheme to solve frac- 
tional Cattaneo equation by a new fractional derivative. Furthermore, they 
analyzed the stability and convergence order of the scheme. The main aim of 
[9] is to prove the existence and uniqueness of the flow of water within a con- 
fined aquifer with Caputo—Fabrizio fractional diffusion for the spatial and the 
temporal variables. In [12], the authors applied the Ritz method with known 
basis functions for a type of Fokker—Planck equation with Caputo—Fabrizio 
fractional derivative. In 2017, Mirza and Vieru [21] proposed the fundamen- 
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tal solutions to time-fractional advection-diffusion equation without singular 
kernel. They applied the Laplace transform and Fourier transforms with re- 
spect to the temporal variable and the space coordinates, respectively. In 
[19], the authors constructed the shifted Legendre polynomials operational 
matrix in order to solve problems with left-sided Caputo—Fabrizio opera- 
tor. A second-order scheme for the space fractional diffusion equation with 
Caputo-Fabrizio is provided in [26]. The main aim of [10] is to solve two 
problems in nonlocal quantum mechanics wherein the nonlocal shrédinger 
equation has been transformed to an ordinary linear differential equation. 
Other interesting papers in the field of the Caputo—Fabrizio derivative are 
found in [1, 2, 3, 5, 8, 11, 18, 14, 20, 22, 23, 25, 17]. 

The main goal of this paper is to derive two new formulas to approxima- 
tion the Caputo—Fabrizio derivative of order a (0 < a < 1). For this pur- 
pose, we use the linear and the quadratic interpolation approximations on 
non-uniform meshes for obtaining the second and the third orders accuracy. 
Besides, we discretize the spatial derivative using the compact finite differ- 
ence scheme. The advantages of the present paper are in the following two 
aspects, that is, the obtained accuracy is independent of the fractional a and 
gives a new high-order accuracy to the time fractional derivative in Caputo-— 
Fabrizio’s sense on non-uniform meshes. In this paper, much attention is 
paid to the numerical aspects. To our knowledge, our interest in Caputo— 
Fabrizio derivative is due to the necessity of using a model describing the 
behaviour of classical viscoelastic materials, thermal media, electromagnetic 
systems, and so on. In fact, the original definition of fractional derivative 
appears to be particularly convenient for those mechanical phenomena, re- 
lated to plasticity, fatigue, damage and with electromagnetic hysteresis. In 
fact, the Caputo—Fabrizio derivative fits to describe material heterogeneities 
and structures with different scales. Hence, we have focused on non-uniform 
meshes. 

The rest of the paper is organized as follows. In Section 2, the derivation 
of the new method on any non-uniform meshes for the Caputo—Fabrizio frac- 
tional derivative of order a (0 < a < 1) in both cases of the second and the 
third order is developed. Three examples are given in Section 3 to support 
the theoretical analysis. Finally concluding remarks are given in Section 4. 


2 Derivation of new method on non-uniform meshes 


2.1 The second-order approximation 


In this section, we focus our attention on deriving the new fractional numer- 
ical differentiation formula in details. By Caputo—Fabrizio fractional deriva- 
tive, we have 
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t 
a a 
fDeu(Dlean = pay fw ex (—ap— Jas 
1 Qf* tn — 
aa [Ww e@en(-o% *)ds, (3) 


where 0 = to < t) < ... < ty = T. We denote the time step by tT, = 
tn — tn-1, 1 <n<QN, and let TMax = MAX1</<N Tl- 

To explain the process of deriving formula, we apply the linear interpola- 
tion polynomial using points (te—1, u(te-1)) and (tx, u(ti.)) for approximat- 
ing u’(s) as follows: 


uth) _ u(tr—-1) 


II,,,u(t) = u(tr) + (t tr) 


th —te-1 
= u(t) + u(t) A () 
and 
TOP (3) aren (5) 


Tk 


Then, the interpolation error formula is given by 


wu" (x) 
2! 


u(t) _ II, ,u(t) = (t _ te-1)(t - tk), tE [tr—1, tr], ERE (th—1, tk), 
1l<k<n. (6) 


Now, substituting (5) into (3), we obtain the L1 formula as follows: 


1 = u(t) = u(te—-1) te tn — Ss 
CF ha = 
0 D; u(t)|r=t,, = Oo 5 _ exp ( = as — ds 


1 
= a (Uk — Up-1) Mz, (7) 
k=1 


where Mp = 4 (cxv( atm—t) — exp ( alerts) ), 


k 
Lemma 1. For any 1 << N, we have M7’ > 0 and M7’, > MP. 


Proof. Note that —“ (tn — tk) > —% (tn — te-1) for 0 < a < 1 and exp(2) is 


l—-a 
a monotone increasing function. Thus one can verify that 


1 th —t tn — tee 
Mg = = (exp a=) exp ( a —))>0 


From the mean value theorem for integrals, we have 
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1 ste th — § tn — & 
Mp=— | exp(- ds = exp (— E (th—1, te): 
k =o ae xp ( os) xp ( aa fy € (th_1, te) 


tn—s 


Clearly, exp ( —-a ns) is a monotone increasing function, then the second 
statement of the lemma follows immediately. 


Theorem 1. Suppose u(t) € C?[0,T]. For any 0 < a < 1, it holds that 


R(ulte))| <a max ee (8) 
Proof. From (2) and (6), we get 
1 = ; — 
R(u(te)) = = [Sf wo — TI «u(s))’ exp (— a 1 “as 
1 n 


I 
—— 
gz 
a 
| 
o 
ua 
7 
a 
w 
nat 
oO 
ia 
Ke) 
— 
| 
Q 
3 
| 
NH 
wn” 
= 


= aoa pe we) (5 t—1)(8 — te) exp ( — at )as 


G@—oy 2 4 


8 Tmax 5 V (to, tn) VRE (tp-1, tk). 


This proves the desired formula. 


We apply the time step over the non-uniform mesh defined as [29] 
™m=(N+1—-n)p, l<n<QN, (9) 
2T 


where L= N(N+I)° 


2.2 The third-order approximation 


Adding an additional point (th-2, u(tk—2)) for k > 2, we obtain a quadratic 
interpolation function T2,,u(t) of u(t) as follows: 
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a 1 u(tr) = u(th—-1) u(trh—-1) — u(tp—2) 
My eu(t) = Th u(t) 4 th — tr-2 | th — tr-1 tr_1 — tr—2 | 
x (t = tk )(t = te—-1) 


= 1 u(tp) = u(th—-1) u(th—-1) = u(th—2) 
- Th wu(t) Tk + Tk—1 | Tk Tk-1 | 
x (t—te)(t—te-1), t © [te-1, te], 


u(t.) — (te 
(Is,,u(t))’ = (ts) _ (fea) | (2t—(te+tx—1)) Aru(te). 
(10) 
For simplicity in what follows, we define: 
He 1 Uk —Uk—-1 — Uk-1 — Uk—2 
Tk + Tk-1 Tk Tk-1 
and 
yl” 
u(t) — Tha e(t) = “OP (ty) (t— tea )(t~ te) 
te [t.—1, trl, Nk © (t._2, th), 2<k<n. (11) 


We substitute (10) into (2) to obtain a new approximation of the Caputo-— 
Fabrizio derivative as follows: 


ty. [ — 
FDO u(t) = | u’(s) exp (-—a— ds 
0 Dy u(t)|t=t, > a (s) p ( ) 


1 a tn — 8 
“ral f (II; ,1u(s)) exp (— a7 — )ds 
n tk 
+> f (Il2,,u(s)) exp (—a——)ds 
kan? th-1 7 
_ il U1 — UO a n s 
=| = [ exp (— a7— )ds 


tk _ bs 
+ 7 eS Arg ae oe) 
tk—1 Tk 
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1 b ae te — 38 
= exp (—a@ ds 
l-a », Tk bed ( ca 
n He tn — 8 
+ So Arte (2s — (th-1 + te)) exp ( ay )ds 
k=2 tea Fa 
1 n 
= Dfu(t)lear,, + D_ BeArur, (12) 
a 
where 
n tn — tk tn —tr-1 
Bi = 2(a-1)| exp( a, ) exp (—a — | 
tn —tey , ty — tr—-1 
tare | exp ( oF) + exp ( ae a )}. (13) 


Moreover, §*D¢ is the L1 method for non-uniform time grid in Section 2, 
where, we define 


Qa Qa 1 . n 
FDFult)lene, =F Deu(t)leata + <5 Y) BeAdte. 4) 
k=2 


wherein, we define a new operator Dfu(t), which is the new fractional nu- 
merical differentiation operator for the Caputo—Fabrizio fractional derivative 
CFT )@ 
5 D*. 


The following lemma states the property of coefficients B?. 


Lemma 2. For any a (0 < a < 1), let 


By = 2(a~1)] exp ( a 


tn — tr 
_— 


+an.|exp(—a 


It holds that 


Bh> BYy_y>-+: > BE > By_4 >--- > By > 0. 


Proof. To prove our statement, we apply the error representation of the trape- 
tn—s 


zoidal formula. For this purpose, we consider the function —2a exp (-a =) 
on the interval [t,-1,t,]. Then 
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th 
Bp = 2a] f exp ( ga 5 (exp ( ao) 
tk—-1 


l-a 


1 th — 8 
= —2a(—75) (exp ( Te l-a es 


2 
a a tn — Ek 
= e tr—-1, tk). 
6G-ap xp (—a = ), && © (te-1, te) 
Since exp (— a=") > 0 and eo > 0, these imply that B? > 0 for 
2<k<n. Besides, exp ( - am") is a monotone increasing function with 


respect to s on [0,7] and this completes the proof. 


Theorem 2. Suppose u(t) € C3[0,7]. For any 0 < a <1, it holds that 


a Ju!" (t)| ju" @ | _s 
; < | ———— ed ; 
ruta) < [Ae mma EON mae MOO re (5) 
Proof. To prove this, we see from (2), (6), and (11) that 
1 " ji th — 8 
Ru(ts)) = =a] (w(s) — Thau(s))' exp (=a =*)as 
1l-al Ji, l-a 
tk 
+ / (u(s) — Il2,nu(s))’ exp (— a——)d 
kao 7 bet 7 
1 tn — $y \t1 
=e (us) — II 1u(s)) exp(—a ina ) . 
——— eee 
=0 
f th —S, @ 
=v (u(s) — Il,,,u(s)) exp (— a = ) 7 ads 
1 . os tn — $y) te 
+ 2L[S [ (uls) — asus)! ew (= 04), 
kao? th-1 
ec 
=0 
” be tn — § Qa 
= ay (u(s) — IIz,,u(s))’ exp (— a i ) ; ds 
han? thoi aa 
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~ (l-a) 2 4 —a 
———$—$_$__ 
TA 
Ww 2 t 
a th — 
Ju LM Te mm) ff exp (—at—)ds ; 
SS 

<ike 


&1 S (to, ti), Vv = (to, tn). 


The second term of last inequality results by the following remark in [29]: 


2 
Remark 1. With regard to max;,<s<t, |(s — tk)(s — te-1)| = is obtained 
at s =tp-1 + 4, consequently we have: 


7 Tk 
(s — ty_2)(s — th_-1) (5 — th) = y (tena a ty—-2) 
re T Te 
= (mat >) S$ lm-a t+ 7%). 


Besides, we know that Tax = maxi<i<n 7. This proves (15). 


3 Numerical applications of the examples 


In the current section, the efficiency of the suggested scheme for the time 
Caputo—Fabrizio fractional diffusion and advection equations are presented 
on three numerical examples in one dimension. The accuracy and the stability 
of the suggested scheme in the paper for different values of M and N are 
tested. In order to carry out our numerical examples, we have used the Maple 
18 software with a PC of 4 GHz CPU and 6 GB memory. The accuracy of 
the proposed scheme is measured by the following error norm 


= . —yNX 
e(N, M) = ey, uz |. 


We denote the numerical convergence orders by 


e(N/2, M) 
Rate =1 —_—_ }. 
ite toes (an) 
Example 1. Suppose 0 < a < 1. Let u(t) = sin(4t). The exact solution is 
obtained from the definition of the Caputo—Fabrizio (2) without the variable 
x. Denote e(N) = |u(tn) — un]. 
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Table 1: Numerical convergence orders in temporal direction for Example 1 


a N_ The second-order The third-order 

e(N) Rate e(N) Rate 

0.5 5 5.6188 x 107? = 2.0699 x 107? - 
10 1.5043 x 107? 1.9012 2.3793 x 1073 ~—- 3.2100 
20 3.9207 x 10-3 ~—- 1.9399 2.6437 x 10-4 —- 3.1699 
40 1.0028 x 10-2 1.9671 3.1297 x 107° = 3.0785 
80 53722. x 10-4 —- 1.9827 3.6530 x 107 = 3.0989 

0.9 5 3.4485 x 1073 — 5.9336 x 107? — 
10 9.3840 x 10-4 1.8777 7.9779 x 10-3 2.8946 
20 2.4994 x 10-* —- 1.9086 1.0192 x 10-3 —- 2.9886 
AO 6.4364 x 107° —- 1.9572 1.2836 x 10-4 2.9892 
80 1.6314 x 107° ~—-1.9801 1.5964 x 107° ~— 3.0073 


From Table 1, this fact is extracted that the computational orders of our 
schemes are independent of the fractional order a. This means that with 
changing values a, the computational orders do not change and the orders 
of convergence of these two cases the second-order and three-order schemes 
should be 2 and 3, respectively. 


Example 2. Consider the time fractional diffusion equation [29] 


Dau(a,t) = Sue) + f(a,t), O<a<1,0<t<1, 
u(x,0) = u°(x), 0<a2<1, (16) 
u(0,t) = B(t), u(1,t)=9(t), O<t<1. 


The exact solution of (16) is u(x,t) = sin(max)t?. 


obtained by substituting u(x,t) into (16). 


The functions can be 


Remark 2. For discretizing the term Dries) 


CFD scheme as follows: 


0 u(a, t) | 62 
Q le=ai = h2 x 
Ox 1 oe 


, we apply the fourth-order 


(xi, t) + O(h*). 


Now, we consider the fourth-order approximation of the second derivative 
of u at point x; as follows: 


07 u(x, t) 62 a 
5) in — h2 Uj; ? 


where u!’ denotes the numerical solution at (x;, tn). Then, a difference scheme 
using third-order formula can be obtained as 
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He Deu? = deur +Hff, 1<i<M-1, 1<n<N, (17) 
ug =P(tn), ut =Plin), ,l<n<QN, (18) 
ul =u (a), O<i< M, (19) 
where 
ate = 7 Oating — Ont g 
Moreover 


Hv; => 


35 (vig + 100; + v4-1), 1<i<M-1, 
Vis 7=O0or M. 


It can be seen that 
h2 
Hv; = (1+ a2) o, 1<i<M-1. 


Lemma 3. (See [31]). consider the function g(x) € C®[ai—1, vi4,], and let 
€(s) = 5(1 — s)? — 3(1—s)°. Then 


9! (@iga) + 10g" (ai) + 9" (@i-a) 
12 
= 4 pl 
oy 70) + 9(@:-1) | [9° (x; — sh) + g° (a; + sh) ]E(s)ds. 


Having seen Tables 2 and 3, we observe that the third-order scheme pro- 
duces better results than the second-order scheme. In the Caputo’s sense, the 
accuracy of the presented method is dependent on a. In this case, the com- 
putational orders for a = 0.4,0.6 and 0.8 are 1.6,1.4 and 1.2, respectively, 
when the theoretical order is 3 — a, whereas the accuracy of the presented 
method is not dependent on the fractional a. Table 4 illustrates the error 
and CPU time of the third-order and the second-order schemes. Having seen 
Table 4, we conclude that the third-order scheme produces more accurate 
results than the second-order scheme. Besides, the third-order scheme needs 
fewer temporal grid size and less CPU time for bigger N. 


Figure 1 exhibits the solution curves at final time T = 1 for different 
values of a = 0.1,0.5 and 0.9 with M = N = 50 for Example 2. Figure 2 
shows the comparison of the absolute errors wherein the third-order scheme 
is more accurate than the second order-scheme. The plots of absolute error 
and the numerical solution for a = 0.1 with M = N = 50 for Example 2 are 
shown in Figure 3. 
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Table 2: Numerical convergence orders in temporal direction with M = 50 for Example 


2 


a N The second-order The third-order 
e(N, M) Rate e(N, M) Rate 
0.25 5 3.6444 x 1074 = 1.9185 x 10-4 = 
10 9.9230 x 107% 1.8768 3.0437 x 1075 = 2.6561 
20. 2.5987 x 1075 1.8299 4.4476 x 1078 = 2.7747 
40 ‘6.6615 x 10~® 1.9639 5.4377 x 1077 — 3.0320 
0.5 5 1.1051 x 10-3 — 5.2234 x 10-4 = 
10 2.9975 x 10-4 1.8823 7.9398 x 10-5 2.7206 
20. -7.8515 x 1075 1.9327 1.1188 x 10-5 2.8272 
40 ~—-2.0126 x 1075 1.9639 1.5315 x 10-& —- 2.8689 
0.75 5 2.6813 x 1073 = 9.2499 x 10-4 S 
10 7.2606 x 10-4 1.8848 1.2632 x 10-4 2.8724 
20. ~=—-1.8975 x 1074 1.9360 1.6398 x 10-5 — 2.9455 
40 = 4.8687 x 10-5 1.9625 2.2262 x 10-8 = 2.8809 


Table 3: Numerical convergence orders in temporal direction with M = 50 for Example 


2 


a N The second-order The third-order 
e(N, M) Rate e(N, M) Rate 
0.4 5 7.4249 x 10-4 = 3.7053 x 1074 7 
10 2.0170 x 10-4 1.8802 5.7480 x 10-5 —- 2.8434 
20. ‘5.2738 x 1075 1.9353 8.0087 x 10-& = 2.8434 
40 ~—-:1.3679 x 107-5 1.9469 1.2630 x 10-7 ~—-. 2.6647 
0.6 5 1.5958 x 107-8 _ 6.9588 x 1074 = 
10. 43238 x 10-* 1.8839 1.0289 x 10-4 2.7577 
20 ~—-1.1306 x 1074 1.9352 1.3961 x 10-5 —-.2.8816 
40 =: 2.8973 x 1075 1.9643 1.7908 x 10-& = 2.9627 
0.8 5 3.1781 x 1073 = 9.4828 x 1074 = 
10 8.6131 x 10-4 1.8836 1.2325 x 10-4. 2.9437 
20 = :2.2530 x 1074 1.9347 1.5675 x 1075 = 2.9751 
40 —-+5.7678 x 10-5 1.9658 1.9745 x 10-& — 2.9889 
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Table 4: The errors and CPU time (seconds) of the third-order and the second-order 
schemes for Example 2 


a The third-order (IM = 50) The second-order (M = 50) 
N e(N, M) CPU(s) N e(N, M) CPU(s) 
0.7 8 9.3694 x 10-4 0.03 5 8.6414 x 10-4 0.09 
24 1.1214 x 10-4 0.06 12 7.2084 x 10-5 0.1 
72 1.2902x 107° = 13.10 24 9.6323 x 107-6 5.75 
0.8 8 1.3174 x 1073 0.06 5 9.4828 x 10-4 0.09 
24 1.5769 x 1074 0.09 10 1.2325 x 1074 0.1 
72 1.8033 107° = 14.30 20 1.5675 x 107° 3.86 
0.9 8 1.9637 x 1073 0.07 4 1.6778 x 1073 0.07 
24 2.3611 x 10-4 0.09 8 1.7168 x 1074 0.09 
72 2.7014x 107° =: 14.53 16 1.9223 x 1075 3.05 


2 03 04 «05 06 07 08 09 1.0 


= Numerical 


Exact 


0+ T T T T T T T T T 
O01 02 03 04 O05 06 07 08 O09 10 


= Numerical 


Exact 


0+ T T T T T 
01 02 03 04 05 


0.6 =40.7 #08 09 1.0 


= Numerical 


Exact 


Figure 1: The solution curves at T = 1 with M = N = 50 for Example 2 
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Figure 2: Comparison of the absolute errors for the second-order (left) and the third- 


order (right) schemes with M = N = 50 for Example 2 


15x 107° 
12x10 

1.x 10° 
75x 10-7 
5.x 1077 
2.5x 10°? 


Numerical sulution 


Absolute error 


Figure 3: Plots of the absolute error (left) with M = N = 50 and the numerical solution 


(right) for Example 2 


Example 3. Consider the time fractional advection equation [29] 


SDeu(a,t) = Me) 4 fie,t), O<e<1,0<t<1, 
u(z,0) = u°(a), 0<a2<1, 
u(0,t) = 8), ul,=e@), 0O<t< 1 


Remark 38. For discretizing the term outs!) 


scheme [29] as follows: 


(20) 


) , we apply the fourth-order CFD 
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Table 5: Numerical convergence orders in temporal direction with a = 0.5 and M = 50 
for Example 3 


N The second-order The third-order 
e(N, M) Rate e(N, M) Rate 
5 1.3158 x 107% - 5.7401 x 1073 — 
10 3.4968 x 1073 1.9118 8.8998 x 10-4 2.6892 
20 9.0919 x 1074 1.9434 1.2614 x 10-4 2.8187 
40 2.3234 x 1074 1.9683 1.7195 x 1076 2.8750 
Ou(a, t) _ Oz 
Ox lien = l a HB 2“ ‘ 


Difference scheme using third-order formula can be obtained as 


AS Deu = dgur + Aff, l1<i<M-1, 1<n<QN, (21) 
Uo — P(tn), uy — p(t), wl < n < N, (22) 
u? = u9 (a4), 0<i<M, (23) 
where 
1 
bg¥i = 5p vit a et), 
and 


Av; = g (dizi + 4; + ¥i-1), 1<i<M-1, 
- Vi; i=Oor M. 


It can be seen that 
h2 
Av; = (1+ 782), 1<i<M-1. 


Like Remark 2 but with slight change, we consider the fourth-order ap- 
proximation of the first derivative of u at point x; as follows: 


Ou(x,t) Og 


gag: = 2 — ila, + OW): 24 
an ln=a; ip Be" )+ O(h*) (24) 


The exact solution of (20) is u(a,t) = sin(rx)t®. Functions can be obtained 
by substituting u(x,t) into (20). 
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The second-order The third-order 


0.85 


0.65 


0.45 


0.254 


T 7 ; ; T : T 0- T T T r 7 r ¥ i T 
03° #04 OS 06 07 O8 09 1.0 0.1 02 03 #04 #O5 06 07 O8 O09 1.0 


= Numerical 


Exact = Numerical Exact 


Figure 4: The solution curves at T = 1 with M = 50, N = 40 for Example 3 


Table 5 confirms that the numerical convergence orders in both the 
second-order and the third-order schemes are close to theoretical results. 
Having seen Table 5, we conclude that the third-order produces better re- 
sults for e(N, M) than that the second-order both in error and accuracy. In 
[29], the accuracy of the presented method is dependent on a. In this case, 
the computational orders for a = 0.5 is 1.5 with the theoretical order of 3—a, 
whereas the accuracy of our schemes is not dependent on the fractional a. 
Figure 4 illustrates the plot of numerical solution and exact solution with 
a =0.5,M = 50 and N = 40 at final time T = 1 for Example 3. 


4 Conclusions 


In the current paper, we have obtained two new fractional numerical differen- 
tiation formulas to approximate the time Caputo—Fabrizio fractional deriva- 
tive of order a (0 < a < 1) on non-uniform meshes. First, the linear and 
quadratic interpolation approximations are considered for the integrand u(t) 
because of obtaining the new formulas. Then, a fourth-order CFD scheme 
is employed for spatial discretization. This difference scheme is led to the 
third-order (second-order) and the fourth-order accuracy in the temporal and 
the spatial variables, respectively. Numerical results are carried out to sup- 
port the convergence orders and show the efficiency of the suggested scheme. 
What distinguishes this paper from our previous studies is its accuracy as- 
pect because the accuracy of the suggested schemes is not dependent on the 
fractional a. 
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